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Abstract 

We compute the next-to-leading order strong interaction corrections 
to gluino-mediated AF = 2 box diagrams in the Minimal Supersymme- 
tric Standard Model. These corrections are given by two loop diagrams 
which we have calculated in three different regularization schemes in 
the mass insertion approximation. We obtain the next-to-leading order 
Wilson coefficients of the AF = 2 effective Hamiltonian relevant for 
neutral meson mixings. We find that the matching scale uncertainty is 
largely reduced at the next-to-leading order, typically from about 10- 
15% to few percent. 
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1 Introduction 

If ongoing and planned experiments discover direct or indirect departures from 
the Standard Model (SM), the next step will be to understand what kind 
of new physics is involved. Detailed predictions for flavour changing neutral 
currents processes in supersymmetry play a crucial role in this program. In 
particular, flavour mixing induced by sfcrmion mass matrices is a pure super- 
symmetric effect with no analogue in the SM and constitutes in general the 
bulk of SUSY contributions to neutral meson mixings. These processes provide 
in turn one of the most sensitive guideline for reconstructing the structure of 
SUSY soft breaking terms. Ultimately, this information will allow us to dis- 
criminate among the many possible mechanisms for SUSY breaking that have 
been proposed in the literature. 

In this paper we present the next-to-leading order (NLO) strong interaction 
corrections to gluino-mediated AF — 2 box diagrams in the Minimal Super- 
symmetric Standard Model (MSSM). We obtain the NLO Wilson coefficients 
of the effective Hamiltonian relevant for neutral meson mixings. We adopt the 
mass insertion approximation [1] which is phenomenologically motivated and 
permits a compact presentation of the results for the Wilson coefficients. 

The complete expressions of the Wilson coefficients at the NLO are col- 
lected in appendix A, where the results are presented in the MS-DRED renor- 
malization scheme. In eqs. (4.7)-(5.2) and (4.8)- (5. 10) we provide the formulae 
required to translate the Wilson coefficients to the MS-NDR and RI-MOM 
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schemes. The relations between the strong couphng constant and the squark 
and gluino masses in the MS-DRED and NDR schemes are given in eq. (A. 3). 

At the LO, strong interaction contributions to AF = 2 processes are de- 
scribed in SUSY by the gluino mediated box diagrams represented in fig. 1. 
These diagrams have been computed in refs. [2]-[3]. NLO corrections to AF = 2 
processes are available for the chargino contributions in the MSSM [4] and for 
the Two Higgs Doublet Model [5] . Both chargino and gluino contributions have 
been then computed in [6] in the MSSM with minimal flavour violation. The 
anomalous dimension matrix for the complete set of four-fermion operators 
entering the effective AF — 2 Hamiltonian has been evaluated at the NLO in 
QCD in refs. [7, 8]. 

This paper completes the NLO determination of the effective Hamiltonian 
by computing the initial conditions for the Wilson coefficients at the high- 
energy supersymmctric scale. Besides the general argument that initial condi- 
tions are needed to obtain scheme-independent results and to achieve NLO scale 
invariance, our calculation is strongly motivated by two additional considera- 
tions. First, the LO coefficients generated by gluino exchange are proportional 
to ttg. Without the NLO computation of matching conditions, it is not pos- 
sible to specify the scale and scheme for the strong coupling, resulting in an 
uncertainty of the LO result much larger than in ordinary weak-interaction 
processes. Second, the new AF = 2 operators generated by gluino exchange 
have surprisingly large anomalous dimensions, so that there is a large scale 
dependence that can only be removed by adding the NLO corrections to the 
matching (see eq. (5.4)). We consider two different regular izat ion schemes 
for ultraviolet (UV) divergences, namely the naive dimensional regularization 
(NDR) and the dimensional reduction (DRED). Infrared (IR) divergences are 
treated both with a gluon mass (in the NDR and DRED schemes) and with 
dimensional regularization (in DRED). The main achievement of the NLO de- 
termination is a strong reduction of the high-energy scale dependence of the 
Wilson coefficients compared to the LO, typically from about 10-15% to few 
percent. Applications of our calculation are studies of B^^g ~ -Bd,s, D — D and 
K — K mixings. Preliminary results for the — B^ mixing case have been 
given in ref . [9] and a complete phenomenological analysis will be presented in 
a forthcoming paper. 

The plan of the paper is the following. In section 2 we introduce the effective 
Hamiltonian approach and the basic formulae used in the matching procedure 
at the NLO. In section 3 we discuss the calculation in the full theory (the 
MSSM) both at the LO and at the NLO. The latter represents the main result 
of the paper. We give details of the calculation and address in particular the 
issues related to the role of evanescent operators in the matching. In section 4 
we present the calculation in the effective theory. The results for the Wilson 
coefficients are discussed in section 5 together with the consistency checks 
between results obtained in the different UV and IR regularization schemes 
and the scaling under the renormalization group equation. Finally, in section 6, 
we draw our conclusions. The complete expressions of the Wilson coefficients, 
both at the LO and at the NLO, are collected in appendix A 
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2 Effective Hamiltonian for AF = 



2 processes 



The effective Hamiltonian for AF = 2 processes in the presence of new physics 
can be written in terms of eight independent four-fermion operators, 

nf/=' = j:Qo, + j2c,d,, (2.1) 

i=l i=l 

where Cj are the Wilson coefficients and we adopt the following basis for the 
local operators Oi 

02 = d'PL V d^PL V , 

03 = d'PL V d^PL b' , 

04 = d'PLb'd^PRb> , 

05 = d'PLb> dPPnh' . (2.2) 

The operators Oi,2,3 ai'c obtained from Oi,2,3 by the exchange L ^ R. The left- 
and right-handed projectors are defined as Pr^l = (l=t75)/2 and 7^ = •j'^Pr^l] 
i,j are colour indices. In eq. (2.2) and in the following we specialized for 
definiteness on the effective Hamiltonian which describes Bd—Bd mixing. In the 
case of Bg, D and K mixings, the replacements {d, b} {s, b}, {d, b} {u, c} 
and {d, b} {d, s} should be considered respectively. 

The evaluation of the coefficients of an effective Hamiltonian involves the 
following two steps: 

1. calculating the amplitude in both the full and the effective theory and 
determining the Wilson coefficients by matching the two amplitudes at 
the high energy scale; 

2. evolving the Wilson coefficients from the high- to the low-energy scale 
where the matrix elements of the local operators can be computed with 
non-perturbative methods, primarily lattice QCD calculations. 

Step 1 depends on the theory under consideration. The new result of this paper 
is the computation of the full theory amplitude in the MSSM up to the NLO 
in the strong interactions. As far as step 2 is concerned, the NLO anomalous 
dimension of the effective Hamiltonian in eq. (2.1) has been calculated in ref. [7] 
and the result confirmed in [8]. 

We now recall the general formulae necessary to perform the matching 
between the full and the effective theories at the NLO. 

The renormalized amplitude in the full theory can be written in the form 

A^uu = E{FP + ^F|''){0^r^ (2-3) 

where (Oi) are the tree level matrix elements of the operators Oi and F^^^ 
and F^^^ represent the LO and NLO contributions respectively. Note that, 
in the case of the AF = 2 SUSY transitions considered in this paper, both 
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F(°) and F^^^ contain an additional factor not factorized out in eq. (2.3). 
It is also worth recalling that, in order to properly normalize the physical 
amplitude, the external quark fields considered to compute the amplitudes 
should be renormalized with their on-shell renormalization constant, defined 
as the pole residue of the quark propagator. In the calculation performed in this 
paper the external fields, as well as the strong coupling constant, renormalize 
differently in the full (MSSM) and in the effective theory, and this gives a 
finite contribution to the matching. In particular, one loop corrections to the 
quark propagator in the full theory include a squark-gluino loop as well as a 
quark-gluon loop, whereas only the latter appears in the low-energy effective 
theory. 

It is convenient to express also the NLO renormalized amplitude in the 
effective theory in terms of tree-level matrix elements of local operators, 

AeS = J:C^{0^)=Y.C^(l + %) {O,)^'^- (2.4) 

By equating the full theory amplitude in eq. (2.3) with the effective one given 
in eq. (2.4) one obtains the expression for the Wilson coefficients at the NLO, 

The functions F^'^^ and r depend in general on the external states. In 
our calculation we have chosen massless external quarks with zero momenta. 
Though this choice considerably simplifies the calculation of the two-loop dia- 
grams in the full theory, it also introduces IR divergences in both the full and 
effective theories, in particular in F^^^ and r. These divergences cancel in the 
Wilson coefficients. Particular care, however, must be taken when regularizing 
IR divergences in dimensional regularization. In this case, the matrix r con- 
tains 1/e poles that give finite contributions to eq. (2.5) once combined with 
both (9(e)-terms entering F^^^ and contributions to F^'^'> of evanescent opera- 
tors. In particular, the summation index k in eq. (2.5) must run in this case 
over both the physical and the evanescent operators, whose specific definition 
will be given in the next section. The evanescent operators which are needed 
instead to define the renormalization scheme of four-fermion operators within 
dimensional regularization are discussed in sec. 4. 

3 Calculation in the full theory 

We now describe the NLO calculation of the Wilson coefficients for AF = 2 
transitions mediated by strong interactions in the MSSM. We will discuss in 
turn the computation of all the elements entering the r.h.s. of eq. (2.5): the 
determination of the LO and NLO amplitudes in the full theory, F^^^ and F^^\ 
is discussed in this section; the calculation of the amplitude in the effective 
theory, expressed by matrix r, will be discussed in section 4. 

As mentioned before, having chosen external quarks with zero masses and 
momenta, the bare amplitudes in both the full and effective theories present 
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UV as well as IR divergences. To regularize both of them we have adopted 
three regularization setups: 

• DRED, with a gluon mass A as IR regulator (DRED-A); 

• DRED, to regularize both UV and IR divergences (DRED-d); 

• NDR, with a gluon mass A (NDR-A). 

The calculation is performed in the mass insertion approximation [1] which is 
phenomenologically motivated and allows a more compact presentation of the 
final results. In order to fix the notation, we recall here the basic formula of 
the mass insertion approximation which provides the expansion of the squark 
mass matrix in the flavour basis around its mean diagonal value, 

(Zt),,(M^),(Z),, = {M%, = Ml[l + ^)i ^ (1 + S)^. . (3.1) 

The matrices Z and M^j are the squark mixing and mass matrix respectively in 
the mass eigenstate basis; M is the squark mass matrix in the super-CKM basis 
{QlQlQlQrQrQr)'i is a mean squark mass, as defined for example in [3]; 
Aij (5ij) are the dimensionful (dimensionless) mass insertions between squarks 
of flavour i and j. We treat Ms as the usual mass parameter in the Lagrangian 
and the 5's as interaction terms. We then expand the AF = 2 amplitude up to 
the second order in the 5's, which provides the first non-vanishing contribution 
in the mass insertion approximation. 

3.1 LO calculation up to 0{e) 

The amphtude of AF = 2 transitions via strong interactions at the LO in the 
MSSM receives contribution from the four box diagrams represented in fig. 1 
for the Bd — Bd mixing case. 

We denote the diagrams in the first and second row of fig. 1 as A-type and 
B-type diagrams respectively and we will extend this notation to the analogous 
topologies entering at the NLO as well (see fig. 2). B-type diagrams entail 
the typical ambiguity in defining the fermion fiow present when dealing with 
Majorana fermions. For a discussion on this point and for the Feynman rules 
of the MSSM wc refer the reader to the refs. [10]- [11]. 

According to eq. (2.5), the Wilson coefficients at the LO are given directly 
by the amphtudes Fj^\ As discussed in the previous section, however, in the 
presence of dimensionally regularized IR divergences, the NLO calculation of 
the Wilson coefficients also requires the evaluation of the LO coefficients of 
the physical operators up to 0{e), as well as the evaluation at the LO of the 
coefficients of the evanescent operators. This is due to the presence of the last 
term in eq. (2.5). In the DRED regularization scheme, one finds the appearance 
of both a d-dimensional metric tensor Qj^i, generated by loop integration (the 
momenta are rf-dimensional) and of a four-dimensional tensor, g^^^, coming 
from the algebra of four-dimensional gamma matrices. Evanescent operators 
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Figure 1: Feynman diagrams describing the gluino contribution to the — 
transition in the MSSM. A cross indicates a mass insertion and the indices h,k,l,m 
label the squark chiralities. The diagrams denoted as A-type and B-type in the text 
are those represented in the first and second row respectively. 



are generated in this scheme by the contraction of Dirac strings with the tensor 
^g^lV1 which can be defined by the following splitting of the metric tensor [12]: 

Qt^v = '^9t,u + i^iiv - ^^M^^ = '^Qt.u + ^Qpiv , (3.2) 

where d = A — 2e and the relations 

9^lu 9p = g,,p , ^g^^,y 9^" = (3.3) 

define the contraction rules in the DRED scheme. The term A^f^,^ is of 0{e) 
and provides our definition of the evanescent operators. In the calculation of 
the LO diagrams we find the appearance of the following evanescent operators: 

= Ag^^d^rLb'd^rR^ , 
^3°^^^° = Ag,,d^rLVd^rRb\ (3.4) 

plus Ef^^^, obtained from i^pRED ^^iq exchange L ^ R. 

We have performed the LO calculation by using the three regularization 
schemes discussed at the beginning of this section. The scheme independent 
results for the Wilson coefficients of the physical operators at the LO, in four 
dimensions, are in agreement with those obtained in ref. [3] and are presented 
for completeness in appendix A. 



3.2 NLO calculation 

The Feynman diagrams entering the calculation of the amplitude at the NLO 
are shown in figs. 2-5. They have been generated by using the Mathematica 
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Figure 2: NLO diagrams generated by gluon corrections to A-type and B-type LO 
topologies. 
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A-type graph 



B-type graph 
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Figure 3: NLO diagrams generated by self-energy corrections to A-type and B-type 
LO topologies. 



A-type graph B-type graph ^ 




Figure 4: NLO diagrams generated by squark corrections to A-type and B-type LO 
topologies. 
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Figure 5: NLO diagrams generated by four squark interaction vertices. 

[13] package FeynArts [14]. The full set of NLO diagrams can be divided in 
four categories. 

1. Gluon corrections, connecting different legs in A- type or B-type LO dia- 
grams. These corrections are collected in fig. 2. 

2. Self-energy corrections of internal legs; these are collected in fig. 3. 

3. Squark corrections, generated by adding to the LO topologies one more 
squark propagator via the quark-squark-gluino interaction vertex. These 
diagrams are shown in fig. 4. 

4. Quartic scalar interactions, generated by the four squark vertex and col- 
lected in fig. 5. 

All these diagrams, except for those belonging to the last category, are gene- 
rated from the LO topologies with the inclusion of an additional loop in all 
possible ways. The last column in figs. 2-5 indicates the number of existing 
diagrams, including the one shown in the figure, that are obtained from the 
latter by performing 90° or 180° rotations around the horizontal, vertical or 
perpendicular axis. Diagrams containing self-energy corrections of the exter- 
nal legs have not been included in the above list. As discussed in section 2, 
however, these corrections have to be taken into account and receive two kinds 
of contributions. QCD contributions mediated by gluons enter the calculation 
of both the full and the effective theory and cancel in the matching, while 
super symmetric squark and gluino corrections give a finite contribution to the 
NLO Wilson coefficients. 

Among the diagrams presented in figs. 2-5, those producing either UV or 
IR divergences are the following ones, 

UV divergent: {A15, A^, ^55^, ^55,, Ajjg, Ajjg, + {A ^ B} 
IR divergent: {Au, A^, Au} + {A ^ B} . (3.5) 

By looking at figs. 2-5 one can see that UV divergent graphs are only those 
containing vertex and self-energy corrections. These graphs provide in par- 
ticular the SUSY contributions to the renormalization of the strong coupling 
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constant and of the squark and gluino fields and masses. IR divergences, in- 
stead, are produced by those diagrams in which a virtual gluon connects two 
external quark lines. These diagrams are in a one-to-one correspondence with 
the diagrams entering the calculation in the efi^ective theory and the whole set 
of IR divergences cancel in the matching. 

We now describe, in some detail, the procedure followed in the evaluation 
of the two-loop diagrams of the full theory. 

Having chosen external quarks with zero masses and momenta, a typical 
two-loop amplitude can be schematically expressed as 



d'^qi d'^q2 FAjqi, q2, fJ',J^, ■ ■ ■) <S) TBiqi, q2, IJ',J^, ■ ■ ■) 
{2ny {2ny {qf - m?)- {q^ - - ^2? - mlf 



where Ta,b represent strings of gamma matrices and loop momenta with sa- 
turated Lorentz indices. To simplify the notation, external quark spinors in 
the amplitude have been omitted. In the denominator, partial fractioning has 
been applied in order to express it in terms of the minimum number of scalar 
propagators, which is equal to three for a two-loop calculation with vanishing 
external momenta. The masses mi^2,3 stand generically for the different masses 
entering the calculation, namely the gluino mass Mg, the mean squark mass 
Mg defined in eq. (3.1) and, when regularizing with a massive gluon, the gluon 
mass A. 

One of the advantages of working with vanishing external momenta is that, 
once the loop integration has been performed, the amplitude in eq. (3.6) turns 
out to be expressed only in terms of strings of gamma matrices, with either 
physical (F^^ ®^^b) or evanescent [E^^ ®E^^) structures, multiphed by scalar 
functions of the particle masses: 



(m) T^2®T^S + H^) E^^®E^S 



(3.7) 



The functions hi{m) are not of interest for our purposes, since the evaluation 
of the Wilson coefficients at the NLO only requires, according to eq. (2.5), 
the projections f/^^ of the two-loop amplitude on the physical operators. The 
complete basis of Lorentz invariant Dirac structures on which we project is 
given by 

^A^B = {7l®7/.l, 7L®7/.fl, Pl®Pl, Pl^Pr, <®<y^^uL}HL ^ R} (3.8) 

where L ^ R indicates the structures obtained by exchanging left and right 

projectors. 

In order to extract directly from a given amplitude V the coefficients 
of the physical operators, we used a basis of orthonormal projectors. These 
are defined as a set of strings of gamma matrices, P^^ Pb\ satisfying the 
orthonormality conditions 



Tr 



r(«) pU) p(i) 



. (3.9) 
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In the DRED scheme the traces are computed in four dimensions. In NDR 
instead, where gamma matrices are d-dimensional objects, the traces are per- 
formed in d dimensions and the orthonormahty conditions (3.9) are required 
to be fulfilled up to and including terms of 0{e); this is sufficient, since the 
two- loop amplitude in the present calculation contains at most 1/e divergences. 
With these requirements the projectors P^^ ® P^^ are uniquely defined. The 
main advantage of using this procedure is that, once the projection is applied 
to an amplitude of the form (3.6), the rcsTilting expression only involves scalar 
integrals. The number of independent two-loop integrations to be performed 
is therefore significantly reduced. 

Besides satisfying eq. (3.9), the projectors must be also orthogonal to the 
evanescent structures. This requirement ensures that, once the projection is ap- 
plied to the r.h.s. of eq. (3.7), no finite contribution coming from the evanescent 
operators is kept in the amplitude. This issue is of relevance in the DRED-rf 
scheme, where IR divergences are dimensionally regularized. In this case, the 
orthogonality of the projectors to the evanescent operators is guaranteed by 
the following observation: all the Dirac structures entering the evanescent op- 
erators in this scheme have uncontracted Lorentz indices and, after the four 
dimensional projections, can only give rise to products of four dimensional Qfj^i^ 
tensors. The latter, in turn, are orthogonal to the evanescent operators in the 
DRED scheme defined as in eq. (3.4), because of the second of eqs. (3.3). 

After the projection has been performed, the evaluation of the two-loop 
integrals is reduced to computing scalar integrals of the form 



This task is greatly simplified by the use of the recurrence relations [15], which 
allow to reduce all scalar integrals of the form (3.10) to a single two-loop master 
integral, I{mi, m^, rris; 1, 1, 1), besides trivial one-loop tadpole integrals.^ The 
result for the master integral T(mi, m2, rris; 1, 1, 1) is given in ref. [19]. 

A further step is required when one of the three masses in the denominator 
of the integral (3.10) is the gluon mass A, introduced to regularize IR diver- 
gences. As a result of having implemented the recurrence relations, one finds 
that the coefficients multiplying the master integral contain negative powers 
of A, up to C(l/A^). The master integral itself must be therefore expanded 
up to (9(A'^). After the expansion, all power divergences must cancel in the 
amplitude and only logarithmic IR divergences remain, which cancel in the 
matching. 

The last step, after the projection and the loop integration, consists in 
expressing the NLO amplitude in terms of tree-level matrix elements of the 
operators in the basis (2.2). This is done by using Fierz rearrangement and 
color algebra. Note, however, that the possibility of expressing the amplitude 
in terms of tree-level matrix elements, up to terms of 0{e), does not occur 

^The application of recurrence relations can be automatically performed by using the 
Tarasov reduction algorithm [16, 17] implemented in the Mathematica program TARCER [18]. 



J(mi, m2, ma; rii, n2, n^) = 



1 



(3.10) 



J {27ry (2^)'^ {ql m;ff^ {ql 



^IT' {{li - <l2f - ml) 
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Figure 6: Feynman diagrams contributing at one loop to the four-fermion operator 
matrix elements in the effective theory. 



diagram by diagram. It only holds, in general, for the complete amplitude. 
This step already provides, therefore, a useful check of the correctness of the 
calculation. 

The sum of the UV renormalized and IR regularized NLO diagrams gives, in 
the notation of eq. (2.5), the functions Fj^\ that represent the main ingredient 
in the NLO evaluation of the Wilson coefficients. 

4 Calculation in the effective theory 

The second step required in the matching procedure is the calculation of the 
amplitude in the effective theory and, in particular, of the matrix r defined 
in eq. (2.4). Using this equation and introducing the renormalization matrix 
Z for the operators Oi, we can write the one-loop matrix elements of the 
renormalized operators as 

We note again that, in the case of the DRED-c? regularization setup, the first 
index i of runs over the evanescent operators too. The reason is that in the 
presence of dimensionally regularized IR divergences the renormalized matrix 
elements of evanescent operators do not vanish. 

Eq. (4.1) shows that the calculation of the matrix r involves two steps: i) 
the determination of the matrix elements of the bare operators (0^)'"^'"^ up to 
one loop and ii) the one loop determination of the renormalization matrix Z. 

As for the calculation of the bare matrix elements, they receive contribu- 
tions in the effective theory only from QCD interactions. The relevant Feynman 
diagrams are those represented in fig. 6, plus the three diagrams obtained by 
performing 180° rotations. Consistency in the matching procedure requires the 
matrix elements in the effective theory to be computed between the same set of 
external states and with the same regularization procedure for IR divergences 
adopted in the full theory. Therefore, we have performed this calculation by 
choosing massless quarks with zero momentum as external states and imple- 
menting separately the three regularization setups: DRED-rf, DRED-A and 
NDR-A. Note, in particular, that the bare amplitudes vanish identically at 
one loop in the DRED-c? scheme, since all loop integrals in this case reduce to 
tadpole massless integrals which vanish in dimensional regularization. 



13 



Eq. (4.1) also indicates that the one loop results for the bare matrix elements 
must be projected onto the basis of the physical operators. This projection 
implies a definition of the evanescent operators. In the DRED rcgularization 
scheme the only evanescent operators entering the calculation are defined to 
be proportional to the tensor Ag^i, of eq. (3.2). Besides the operators specified 
in eq. (3.4), we also find the appearance of the evanescent operators 

= A^^, dV^ft^- (4.2) 

In the NDR scheme, instead, both Dirac and Fierz evanescent operators 
must be introduced. Dirac evanescent operators arc defined from the orthogo- 
nality condition to the Dirac projectors (see eq. (3.9)), 

= . (4.3) 

The complete list is given in rcf. [8]. As for the Fierz evanescent operators, 
they are defined without introducing in the four dimensional Fierz relations 
arbitrary terms of 0{s); for example, the 7y^L (8) 7^ Fierz evanescent operator 
reads 

E^"" = d%LV d^rL^' - d%LV d^itV (4.4) 

and similarly for the other gamma structures. 

According to eq. (4.1), the second ingredient in the determination of the 
matrix r is the one-loop calculation of the renormalization matrix Z. This 
requires again the evaluation of the Feynman diagrams shown in fig. 6. In 
this case however, in order to identify the UV divergences within dimensional 
rcgularization, one can either regularize the IR divergences with a fictitious 
gluon mass or consider a set of IR finite external states, for instance off-shell 
quarks with fixed momentum p. 

In both the MS-DRED and -NDR rcgularization schemes the renormaliza- 
tion matrix of the physical operators is determined by applying the modified 
minimal subtraction prescription. Evanescent operators, instead, must satisfy 
a different renormalization condition. For IR finite configurations of external 
states, this condition reads 

(Ei(/i)) =0 in the limit ci ^ 4 , (4.5) 

and holds at any value of the renormalization scale n [20]. It guarantees that 
the evanescent operators do not play any role when going back to four dimen- 
sions and can be eventually removed from the operator basis of the effective 
Hamiltonian. 

The final result for the matrix r defined in eq. (4.1) depends on several 
choices done in the calculation: the external states, the IR regulator (when IR 
divergences are present) and the renormalization scheme of the local operators. 
Thus we end up with three different matrices, ^dred-a ^^^^^ ^ndr-a_ 

Here we only present the results for the differences Ar between these matrices 
because, at variance with the r's, they are independent of the specific choice of 



Tr 



7T.(«) p(i) 7T.W pli) 
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both the external states and the IR regulator. As can be seen from eq. (4.1), the 
matrices Ar provide the relation between operators renormalized in different 
schemes. In the case of the MS-NDR and DRED schemes, for instance, this 
relation reads 



^MS-NDR _ / ]^ _|_ A^^NDR/DRED \ ^^ ^MS-DRED g'^ 

where Ar'^^^/^^^^ = r^^^ - r°^^°. For this matrix we obtain the result 



^^NDR/DRED _ 



/ -3 











\ 





-13/3 


-1/3 











-29/6 


7/6 

















-5/3 


-3 


V 








-7/2 


-1/6 y 



(4.7) 



in the basis Oi, . . . , C5 of eq. (2.2). Since chirality is conserved by QCD in- 
teractions in the limit of massless quarks, the corresponding matrix for the 
operators Oi,2,3 is equal to the 3x3 submatrix for Oi,2,3 in eq. (4.7) and the 
two sets of operators do not mix. 

In addition, wc provide the matrix Ar relating the MS-DRED with the so 
called RI-MOM scheme in the Landau gauge [21]. This is useful because this 
scheme is frequently used in lattice QCD calculations of the hadronic matrix 
elements. This matrix reads: 



^^DRED/RI ^ 



/ -| + 8 ln2 





\ 

^2x2 
-62x2 J 



(4.8) 



with 



and 



12x2 



67 
9 
28 



9 9 



f ln2 
28 ln2 



-l + f ln2 
-f + f ln2 



2x2 



13 



^ + 2 1n2 



ln2 1 + 2 ln2 



-|-|ln2 



(4.9) 



(4.10) 



The results in eqs. (4.7) and (4.8) can be also combined to obtain the matrix 
relating the MS-NDR with the RI-MOM scheme: Ar^^^/^^ = Ar^°^/°^^° + 

^^DRED/RI 



5 Results and checks of the calculation 

In the previous sections we have described the calculation of the two ingredi- 
ents needed to obtain the Wilson coefficients at the NLO: the LO and NLO 
amplitudes in the full theory, F^^^ and F^^\ and the matrix r in the effective 
theory. The NLO Wilson coefficients are finally determined using eq. (2.5). 
They bear a dependence on both the renormalization scheme and scale. These 
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dependences only arise at the NLO and allow one to perform useful checks 
of the calculation. The relations among the results for the coefficients as ob- 
tained in the three regularization setups, DRED-A, NDR-A and DRED-rf, will 
be discussed in the following subsection. The scale dependence of the Wilson 
coefficients must satisfy the renormalization group equation, and this constraint 
will be addressed in subsection 5.2. 

5.1 Regularization and renormalization scheme depen- 
dence 

The results for the coefficients obtained in the DRED-A setup must be equal 
to those obtained in DRED-d, since the Wilson coefficients cannot depend on 
the IR regulator. Indeed, upon explicit comparison, they are found to be in 
agreement. We emphasize that this is a non-trivial check of the calculation. 
Indeed, whereas the computation in the DRED-A scheme presents basically no 
subtlety, the one in the DRED-d regularization entails the inclusion in the full 
theory of the LO contributions up to 0{e) and of the evanescent operators. 
All these contributions should sum up to reconstruct the results obtained by 
using the gluon mass as IR regulator. 

The results for the Wilson coefficients obtained in the MS-DRED and NDR 
renormalization schemes differ because the coefficients arc scheme dependent 
quantities. They can be compared using the scheme independence of the effec- 
tive Hamiltonian: 

{b,d\n,s\b,d) = ^Cf^^°(6,rf|g,|6,5)°^™ = Cf^^(6,rf|g,|6,d)^'^^ . (5.1) 

i 

The relation between renormalized operators in the two schemes has been writ- 
ten in eq. (4.6) in term of the matrix Ar. From eq. (5.1), it then follows that 
the same matrix also relates the coefficients in different schemes, e.g. 

Cf °^(M„ M-„ a,) = E (l - ^Ar^°^/™) C™(M„ M-„ a,) , (5.2) 

where /^t-NDR/dred given in eq. (4.7). Notice that in eq. (5.2) the transposed 
matrix Ar'^ enters. The coupling constant and the SUSY masses Mg and 
Mg in the previous equation are also scheme dependent quantities. This de- 
pendence starts at 0{as) and must be taken into account in the matching at 
the NLO. In order to verify eq. (5.2), therefore, one needs to express all the 
couphngs in the same scheme. The required relations are [22]: 

^NDR^^DRED (^^ + ^ _ c<^)^ 

Mj^^^ = MP^^d (l + ^7V^) (5.3) 

Mf°^ = M™(l + 0(a^)) , 

where = 3 and Cp = 4/3 are the SU(3)c color factors. The strong coupling 
constant cxg in eq. (5.3) indicates the coupling of the quark-squark-gluino ver- 
tex. A different relation is found for the quark-quark-gluon coupling, which 
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differs from dig in the NDR scheme because this regularization breaks supersym- 
metry [22] (sec cq. (A. 3)). In the present calculation, since the shifts expressed 
by eq. (5.3) arc 0{as), they have to be implemented only in the LO amplitude, 
where only the coupling appears. We find that our results for the Wilson 
coefficients as obtained in the DRED and NDR schemes consistently satisfy 
eq. (5.2), with the matrix Ar^°^/°^^° given in eq. (4.7). 



5.2 Renormalization scale dependence 



Beyond LO, the Wilson coefficients acquire an explicit dependence on the renor- 
malization scale /I. This dependence is controlled by the renormalization group 
equation, which provides therefore an additional check of the calculation. 

The renormalization group equation for the Wilson coefficients of the MSSM 
[23, 24] can be written as 



d 



+ 



da. d 



+ 



dM? 



d 



9 In /x^ d In //2 dag d In jj? dM? 



dM^ d ^ dAx d 



C{ix) = , 



(5.4) 



and takes into accoimt the scale dependence of all the quantities entering the 
coefficients, namely the strong coupling constant a^, the squark and gluino 
masses Mg and Mg and the dimensionful mass insertions A^, with X — 
LL^ RRj LR , RL. 

The matrix 7 in eq. (5.4) is the anomalous dimension matrix of the four- 
fermion operators (2.2) in the effective theory {i.e. QCD). It can be expanded 
as 

7(a.) = Z-^3^ = ^7o + 0(«D (5.5) 



din// 



where, for the LO anomalous dimension 70, we obtain the expression 

\ 



7o = 



/ 4 









-28/3 
16/3 









4/3 
32/3 











-16 

-6 2/ 



(5.6) 



in agreement with eqs. (11)- (13) of ref. [25]. 

The renormalization group equation for the strong coupling constant in the 
MSSM reads 



MSSM, 



da<. 



dln/i^ 



_g^MSSM^^(^3)^ 



(5.7) 



with /^o^SSM _ 3jy^ _ jy^^ 

The scale dependence of the squark and gluino masses, Mg and Mg, is 
described instead by the equations 



1 dMf 
Mfdln^j? 



as 

4:71 



(5.8) 



17 



(^) |Ci" / Ci 



4t 

5 

3 

5 

2 

5 






■ r 







P P LR p LL p LR p LL 

'^4 '-'4 



Figure 7: Comparison between the LO and NLO contributions to the MS-DRED 
Wilson coefRcients at the scale /x = Ms and at the reference value Mg = Ms- For 
the coefficients C4 and C5 the contributions proportional to Sll^rr and SlrSrl are 
shown separately. 



where 7]^"; = ACpMl/M^^ and 7!"! = 2(3f^'^'^ = 2(3Ar, - A^^). 

Finally, the running of the dimensionful mass insertions Ax is expressed by 

O + OK^) , 

-^if ^LRiRL) + 0{al) (5.9) 
with 7X"^ = 2Cf. 

By using the results given in eqs. (5.7)-(5.9), we have then verified that our 
expressions for the Wilson coefficients exhibit at the NLO the correct renor- 
malization scale dependence predicted by eq. (5.4). 

5.3 Discussion of the results 

We conclude this section by presenting and discussing the final results ob- 
tained for the Wilson coefficients at the NLO. The complete expressions of 
these coefficients, in the MS-DRED renormalization scheme, are collected in 
appendix A. 

In order to illustrate the typical size of the computed NLO corrections, we 
show in fig. 7 the values of the NLO contributions to the Wilson coefficients 
normalized to their expected size, namely the corresponding LO coefficients 
multiplied by as{Ms)/T^. For the purpose of illustration, in this comparison 
we set the scale = Mg and put Mg = Mg. As can be seen from the plot, 
in several cases the NLO coefficients turn out to be larger than what naively 
expected. Of course, this conclusion applies to the MS-DRED coefficients and 
could change in a different renormalization scheme. 
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The Wilson coefficients depend on tfie matdiing scale which can be cho- 
sen around a typical SUSY scale, e.g. the average squark mass M,. An im- 
portant achievement of the NLO calculation is a significant reduction of this 
dependence with respect to the LO approximation. This is illustrated in fig. 8 
where we show the LO and NLO predictions for the i^d-niesons mass differ- 
ence Am^ as a function of the high-energy scale ^l chosen for the matching. 
These predictions are obtained by adding the SUSY contributions to the ref- 
erence SM value, Am^^ = 0.496 ps~^. The hadronic matrix elements are 
evaluated by using the lattice QCD results of ref. [26] for the B-parameters 
and fs^ = 189 MeV. We set Ms = Mg = 350 GeV and consider two cases for 
mass insertion coefficients, 5ll — Srr — 0.12 exp[— i 23°] (upper plot in fig. 8) 
and Sir — Srl — 0.026 exp[—i 23°] (lower plot in fig. 8), chosen to give a SUSY 
contribution compatible with the present measurement taking into account the 
SM uncertainty. Clearly, the reduction of the scale dependence found at the 
NLO quantitatively depends on the specific values chosen for the mass insertion 
parameters. 

Prom fig. 8 we see that the SUSY prediction for Am^ varies, at the LO, by 
approximately ±16% (±8%) in the LL/RR (LR/RL) case, when the scale fi is 
varied in the typical range between Ms/2 and 2Ms. With the NLO calculation, 
the dependence on the matching scale is reduced by a factor two or more, i.e. 
at the level of ±5% (±2%) percent. 

We conclude this section by observing that phenomenological applications 
require the knowledge of the hadronic matrix elements. These are usually 
computed on the lattice where, in order to perform a fully non-perturbative 
renormalization, the RI-MOM renormalization scheme [21] is needed. This is 
the scheme adopted for instance in refs. [27, 28] and [26], where lattice results 
for the complete basis of four-fermion operator matrix elements relevant for 
K — K and B^^g — Bd,s systems have been presented. In these cases, the results 
for the Wilson coefficients given in appendix A in the MS-DRED scheme must 
be converted to the Rl-MOM scheme. This can be easily done using the relation 
analogous to eq. (5.2), namely 



The matrix (Ar^^^'^/^^)^ which performs the matching between the two schemes 
can be obtained by transposing the matrix in eq. (4.8). 

6 Conclusions 

In this work we have computed the NLO strong interaction corrections to the 
Wilson coefficients relevant for AF — 2 transitions in the MSSM with the 
mass insertion approximation. The complete expressions for the coefficients 
are given in appendix A in the MS-DRED scheme. We also give in eqs. (5.2) 
and (5.10) the formulae required to translate the Wilson coefficients at the 
NLO from the DRED to the NDR and the Rl-MOM renormalization schemes, 
which might be useful for phenomenological applications. 




(5.10) 
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Figure 8: LO and NLO predictions for the Bd-mesons mass difference Am^ obtained 
by adding the SUSY contributions proportional to SlLi^rr (top) and 5lr, Srl (bot- 
tom) to the SM prediction. See text for the reference values of the input parameters. 
The results are plotted as functions of the matching scale /i. 
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Theoretically, the NLO calculation of the Wilson coefficients is required to 
cancel the corresponding renormalization scale and scheme dependence of the 
renormalized operators. Once combined with the NLO anomalous dimension 
of the four-fermion operators given in ref. [7], our results allow to perform a 
complete NLO analysis of AF = 2 transitions in the MSSM. The phenomeno- 
logical analysis will be presented in a forthcoming publication. In this study we 
have shown that, by considering as a reference example the theoretical predic- 
tion of the Bd-raeson mass difference Arrid, the uncertainty due to the choice 
of the the high-energy matching scale is largely reduced going from the LO to 
the NLO, typically from about 10-15% to few percent. 
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A Results for the Wilson coefficients 

In this appendix we collect the complete NLO expressions of the Wilson coef- 
ficients entering the effective Hamiltonian which describes AF = 2 transitions 
mediated by strong interactions in the MSSM. 

Wc consider the complete basis of four fcrmion operators given in eq. (2.2). 
The coefficients for the operators (5i,2,3 are obtained from those of the operators 
Qi,2,3 by simply exchanging L <-> it! in the mass insertion parameters. For this 
reason, we will not present in the following their explicit expressions. 

The Wilson coefficients are written as 



Q(/.) = cf)(/x)+cr^(^) 



(1), 



(A.l) 



where fi is the scale used in the matching procedure. To simplify the notation, 
in the following, we will not write explicitly the /i dependence of masses and 
couplings. 

The LO calculation of the Wilson coefficients has been performed in ref. [3] 
and we agree with their results. These coefficients read 
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/5 107a; lla;^ 19x^ /11a; 13a;2\ , \, , 1 

where x = M?/M^ with Mg the gluino mass and Mg the average squark mass. 
The dimensionless mass insertion parameters are understood to be 5^2, ^13, ^23 
and for the cases of K, Bd, Bg and D mixings respectively. 

At the NLO, the Wilson coefficients are scheme dependent quantities. Here 
we present the results for the operators renormalized in the MS-DRED scheme, 
and we adopt the same scheme also for the strong coupling constant and 
for the squark and gluino masses. Eq. (5.10) can be then used to convert 
the coefficients to the RI-MOM scheme, frequently adopted in lattice QCD 
calculations of the corresponding hadronic matrix elements. As for the strong 
coupling constant and the squark and gluino masses, they can be converted to 
their counterparts in the MS-NDR scheme by using [29, 22] 



^NDR ^ ^DRED , ^ 



Att 3 



= mP^^° (1 + ^N,^ (A.3) 

^NDR ^ ^DRED ^ ^^^2 



We now present the NLO expressions for the Wilson coefficients. In the 
following, the symbol Li2(x) denotes the dilogarithm function defined as 

Li2(x) = -/ ^ln(l-t) (A.4) 
Jo t 

and all couplings and masses are understood to be renormalized at the same 
scale 11. The coefficients read 
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